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Summary 

Unlike  typical  earthquakes,  rockbursts  and  mine  collapses  often  have  regional  waveforms  that  may  be 
confused  with  nuclear  explosions.  Such  seismic  signals  could  present  false  alarms  to  organizations  tasked 
with  verifying  a  Comprehensive  Test  Ban  Treaty.  Mines  may  also  provide  potential  locations  for 
evasion,  through  partial  decoupling  or  triggered  shear  failure.  The  Coeur  d’Alene  district  in  northern 
Idaho  is  the  most  rockburst-prone  mining  district  in  the  United  States,  producing  a  daily  average  of  about 
3  events  with  magnitudes  ranging  from  0  to  4.  The  hypocenter  locations  of  many  of  these  events  are 
determined  to  an  accuracy  of  within  30  meters  by  in-mine  rockburst  monitoring  systems  operated  by  the 
»  mining  companies.  The  largest  and  most  frequent  events  are  located  near  an  area  of  active  mining  at  a 

depth  of  1650  meters  at  one  of  three  producing  mines.  Focal  mechanism  studies  based  on  these  systems, 
and  a  district-wide  surface  seismic  network,  indicate  dominantly  dilatational  first  motion  events.  We 
deployed  a  local  seismic  system  in  the  summer  of  1998  to  determine  the  source  parameters  and  moment 
tensors  for  seismic  events  in  the  region,  in  particular  to  confirm  the  shear  implosional  mechanisms 
suggested  by  the  first-motion  studies. 

Characterization  of  mining-related  seismic  events  is  critical  to  the  implementation  of  the  Comprehensive 
Nuclear  Test-Ban  Treaty  monitoring  system.  Mine  seismicity  from  uncontrolled  sources  are  particularly 
troublesome,  with  magnitudes  up  to  4  or  5,  producing  signal  strengths  comparable  to  1-10  kiloton 
contained  nuclear  explosions.  In  this  study,  ground  motions  from  four  surface  stations  and  one 
underground  station  were  inverted  to  determine  the  complete  moment  tensor  of  a  relatively  large  rock 
burst  in  the  Coeur  d’Alene  Mining  District  in  northern  Idaho.  Many  rock  bursts,  induced  by  the  mining 
process,  show  shear  implosional  source  mechanisms,  a  feature  that  should  be  useful  to  discriminate 
mining-related  events  from  underground  explosions. 
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kilo  pascal  (kPa) 

bar 

1.000  000  x  E  +2 

kilo  pascal  (kPa) 

bam 

1.000  000  X  E  -28 

meter2  (m2) 

British  thermal  unit  (thernochemical) 

1.054  350  x  E  +3 

joule  (J) 

calorie  (thermochemical) 

4.184  000 

joule  (J) 

cal  (thermochesriical/cnt2) 

4.184  000  x  E  -2 

mega  jcule/m2  (MJ/trt2) 

curie 

3.700  000  X  E  +1 

*giga  bacquerel  (GBq) 

degree  (angle) 

1.745  329  x  E  -2 

radian  (rad) 

degree  Fahrenheit 

tk  -  (t°f  +  459. 67) /l. 8 

degree  kelvin  (K) 

electron  volt 

1.602  19  x  E  -19 

joule  (J) 

erg 

1.000  000  x  E  -7 

jcule  (J) 

erg/seoond 

1.000  000  x  E  -7 

watt  (W) 

foot 

3.048  000  x  E  -1 

meter  (m) 

foot-pcund-force 

1.355  818 

jcule  (J) 

gallon  (U.S.  liquid) 

3.785  412  x  E  -3 

meter3  (m3) 

inch 

2.540  000  x  E  -2 

meter  (m) 

jerk 

jcule/kilogram  (J/kg)  radiation  dose 

1.000  000  x  E  +9 

jcule  (J) 

absorbed 

1.000  000 

Gray  (Gy) 

kilotcns 

4.183 

tera  joules 

kip  (1000  lbf) 

4.448  222  x  E  +3 

newton  (N) 

kip/inch2  (ksi) 

6.894  757  x  E  +3 

kilo  pascal  (kPa) 

ktap 

1.000  000  X  E  +2 

newtcn-seccnd/m2  (N-s/m2) 

tnicrcn 

1.000  000  X  E  -6 

meter  (m) 

mil 

2.540  000  x  E  -5 

meter  (m) 

mile  (international) 

1.609  344  x  E  +3 

meter  (m) 

ounce 

2.834  952  x  E  -2 

kilogram  (kg) 

pound- force  (lbs  avoirdupois) 

4.448  222 

newton  (N) 

pound-force  inch 

1.129  848  x  E  -1 

newton-meter  (N-m) 

pound- force/ inch 

1.751  268  x  E  +2 

newton/meter  (N/m) 

pound-  force/foot2 

4.780  026  x  E  -2 

kilo  pascal  (kPa) 

pound-  force/inch2  (psi) 

6.894  757 

kilo  pascal  (kPa) 

pound-mass  (Ubm  avoirdupois) 

4.535  924  xE-1 

kilogram  (kg) 

pcund-mass-foot2  (moment  of  inertia) 

4.214  011  x  E  -2 

kilogram-meter2  (kg-m2) 

pcund-mass/foot3 

1.601  846  x  E  +1 

kilogram-meter3  (kg/m3) 

rad  (radiation  dose  absorbed) 

1.000  000  X  E  -2 

♦♦Gray  (Gy) 

roentgen 

2.579  760  x  E  -4 

cculocrb/kilogram  (C/kg) 

shake 

1.000  000  x  E  -8 

second  (s) 

slug 

1.459  390  X  E  +1 

kilogram  (kg) 

torr  (ttm  Hg,  0°  C) 

1.333  22  x  E  -1 

kilo  pascal  (kPa) 

*  *The  bacquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  *  1  event/s. 

**The  Gray  (GY)  is  the  SI  unit  of  absorbed  radiation. 
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SECTION  1 


DEPLOYMENT  OF  A  LOCAL  NETWORK  AND  SEISMIC  DATA  COLLECTION  OF 

NORTHERN  IDAHO  ROCKBURSTS 


1.1  Summary  Statement. 

Unlike  typical  earthquakes,  rockbursts  and  mine  collapses  often  have  regional  waveforms  that  may  be 
confused  with  nuclear  explosions.  Such  seismic  signals  could  present  false  alarms  to  organizations  tasked 
with  verifying  a  Comprehensive  Test  Ban  Treaty.  Mines  may  also  provide  potential  locations  for 
evasion,  through  partial  decoupling  or  triggered  shear  failure.  The  Coeur  d’Alene  district  in  northern 
v,  Idaho  is  the  most  rockburst-prone  mining  district  in  the  United  States,  producing  a  daily  average  of  about 

3  events  with  magnitudes  ranging  from  0  to  4.  The  hypocenter  locations  of  many  of  these  events  are 
determined  to  an  accuracy  of  within  30  meters  by  in-mine  rockburst  monitoring  systems  operated  by  the 
mining  companies.  The  largest  and  most  frequent  events  are  located  near  an  area  of  active  mining  at  a 
depth  of  1650  meters  at  one  of  three  producing  mines.  Focal  mechanism  studies  based  on  these  systems, 
and  a  district-wide  surface  seismic  network,  indicate  dominantly  dilatational  first  motion  events.  We 
deployed  a  local  seismic  system  in  the  summer  of  1998  to  determine  the  source  parameters  and  moment 
tensors  for  seismic  events  in  the  region,  in  particular  to  confirm  the  shear  implosional  mechanisms 
suggested  by  the  first-motion  studies. 


1.2  Introduction. 

Rockbursts  are  seismic  events  that  cause  visible  damage  to  underground  excavations.  Rockbursts  are 
major  hazards  in  mines  worldwide,  particularly  in  Canada  and  South  Africa,  where  they  cause  production 
delays,  losses  of  millions  of  dollars  per  year,  in  addition  to  injuries  and  death  to  miners.  Additionally, 
Heuze  (1994)  pointed  out  the  potential  to  use  these  events  to  camouflage  nuclear  tests.  Considerably  more 
study  is  needed  to  characterize  systematic  differences  between  mine  tremors  and  other  seismic  sources. 
Mine  explosions  and  rockbursts  are  similar  in  some  ways  to  natural  earthquakes  in  that  they  are  often 
Telated  to  shear  failures  on  fault  planes  in  a  rockmass  (McGarr,  1984;  McGarr,  1971;  Spottiswoode  and 
McGanr,  1975;  McGarr  et  al.,  1979).  However,  results  from  focal  mechanism  studies  of  mining  events 
indicate  that  alternative  mechanisms  other  than  simple  shear  failure  are  common. 

Simple  explanations  exist  for  the  deviation  of  the  moment  tensor  solutions  from  double-couple  solutions, 
including  curved  failure  planes,  incorrect  hypocenter  depths,  and  multiple  ruptures.  However,  in  the  case 
of  mining  events,  a  volumetric  component  seems  to  be  common,  as  is  indicated  by  a  preponderance  of 
■dilatational  P  wave  arrivals  (Williams  and  Arabasz,  1989;  Wong  et  al.,  1989;  Stickney  and  Sprenke, 
1993).  Shear  implosional  models  for  mining  events  have  been  proposed  by  Teisseyre  (1980),  and  Rudejej 
and  Sileny  (1985).  At  the  other  extreme,  shear-tensional  models  for  mine  events  have  also  been  proposed 
(Teisseyre,  1980;  Gibowicz  and  Kijko,  1994). 

The  failure  mechanisms  associated  with  mining-related  seismic  events  in  the  Coeur  d’Alene  mining 
district  of  northern  Idaho  are  uncertain.  Problems  associated  with  rockbursting  cause  accidents  and 
increase  the  cost  of  mining.  The  first  step  toward  mitigating  rockburst  hazards  is  a  fundamental 
understanding  of  their  cause.  Prior  to  1980,  rockbursts  in  the  district  were  attributed  to  the  failure  of 
pillars  ahead  of  advancing  mine  faces;  during  the  1980s  however,  it  became  apparent  from  observations 
of  in-  mine  damage  that  shear  slip  along  faults  and  bedding  planes  was  also  a  source  of  rockbursts  (Blake, 
1984;  Blake  and  Cuvelier,  1990).  Several  rigorous  studies  of  mine  tremors  in  deep  South  African  gold 
mines  suggest  that  natural  earthquakes  and  mining-  related  seismic  events  share  a  similar  focal 
mechanism:  pure  shear  failure  (McGarr,  1971, 1984;  McGarr  et  al.,  1975;  Spottiswoode  and  McGarr, 
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Figure  1.  Location  Map  of  Mines  in  the  Coeur  d’Alene  Mining  District  of  Northern  Idaho. 

1975).  However,  the  presence  of  void  space  near  mine-related  foci  theoretically  also  allows  for  a 
compressive  "implosional"  component  of  stress  release  (Teisseyre,  1980;  Stiller  et  al.,  1983;  Kozak  and 
Sileny,  1985;  Rudajev  and  Sileny,  1985).  Many  researchers  have  reported  a  preponderance  of  dilatational 
first  motions  associated  with  mining-  related  seismic  events:  in  South  African  gold  mines  (Gane  et  al., 
1952;  etc.),  and  in  coal  mines  worldwide  (See  Wong  and  McGarr  (1990)  and  Wong  et  al.  (1989)  for 
reviews).  However,  in  a  review  of  the  focal  sphere  coverage  of  these  observations,  Wong  and  McGarr 
(1990)  could  find  only  five  cases  of  implosional  failure  that  were  worthy  of  consideration  and  they 
discredited  four  of  these  five  cases  on  the  bases  of  inadequate  focal  sphere  coverage  or  inaccurate 
hypocenter  locations.  They  also  found  insufficient  evidence  to  validate  a  shear-implosional  model  for 
coal  mine  slip  events. 

Rockbursts  and  mine  collapses  are  often  observed  to  have  regional  waveforms  that  are  unlike  typical 
earthquakes  and  which  can  be  confused  with  nuclear  explosions  (Taylor,  1994,  Bennett  et  al.,  1995, 
Bennett  et  al.,  1996).  Mining-induced  events  cannot  be  discriminated  on  the  basis  of  hypocentral  depth. 
The  relatively  low  excitation  of  Love  and  Rayleigh  waves  tend  to  make  mining-induced  events,  especially 
those  with  a  large  implosional  component  to  the  moment  tensor,  appear  similar  to  explosions  when  using 
classical  Ms  vs.  mb  discriminants.  Bennett  et  al.,  (1995,  1996)  indicate  that  the  S/P  or  Lg/P  ratios  for 
mining-induced  events  around  the  globe  are  larger  than  observed  for  nuclear  explosions, which  makes  the 
regional  signals  appear  more  similar  to  typical  earthquakes.  Conversely,  Taylor  (1994)  observed  a 
mining-related  event  in  Utah  that  had  a  low  Lg/Pg  ratio,  and  other  signal  characteristics  including  a  lack 
of  Love  waves,  that  caused  this  event  to  be  classified  as  an  explosion  in  an  extensive  discrimination  study 
in  the  western  U.S.  (Taylor,  1991).  Pechman  et  al.  (1995)  note  that  a  mine  collapse  in  southeastern 
Wyoming  also  had  mb-Mw  and  ML-Mw  differences  that  were  more  typical  of  nuclear  explosions,  and 
also  noted  the  lack  of  Pnl  and  Love  waves,  also  noted  to  be  a  feature  of  nuclear  explosion  waveforms. 
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13  Previous  Work  in  the  Coeur  d'Alene  District. 


The  Coeur  d'Alene  mining  district,  one  of  the  largest  lead-,  zinc-,  silver-producing  areas  of  the  world,  lies 
along  the  Osbum  fault  in  northern  Idaho  near  the  Montana  border  (Figure  1).  The  tectonic  setting  of  the 
district  within  the  larger  Lewis  and  Clark  Zone  has  been  described  by  Hobbs  et  al.  (1965);  Wallace  et  al. 
(1990);  and  Stickney  and  Bartholomew  (1987).  Although  seismicity  in  the  district  has  historically  been 
related  to  mining  activity,  tectonic  events  also  occur  (Sprenke  et  al.,  1991;  Sprenke  and  Stickney,  1993; 
Breckenridge,  Sprenke,  and  Stickney,  1995).  To  gather  more  information  on  the  source  mechanisms  of 
mining-related  seismic  events  in  the  Coeur  d'Alene  district,  a  district-wide  seismograph  network  was 
operated  using  portable  instruments  for  one  week  during  September  1990  (Stickney  and  Sprenke,  1993), 
and  a  district-wide  network  using  up  to  28  stations  from  1992-1994  (Jung  et  al.,  1994).  As  a  result  of 
these  studies,  we  have  observational  evidence  that  dilatational  first  motions  are  prevalent  in  this  deep 
metal-mining  district,  probably  the  result  of  implosional  and  shear-implosional  source  mechanisms 
involving  mine  openings. 

Using  data  from  the  district-wide  networks,  hypocenter  locations  determined  by  the  in-mine  systems, 
were  used  to  construct  a  travel  time  curve  (Figure  2)  for  the  district  using  impulsive  body  wave  travel 
times  for  19  Lucky  Friday  Mine  and  two  Galena  Mine  events.  The  results  indicate  a  half  space  velocity 
model  with  a  P  wave  velocity  of  5.74+/-0.02  km/s  and  a  uniform  S  wave  velocity  of  3.47+/-0.05  km/s 
(Stickney  and  Sprenke,  1993). 

The  preferred  modes  of  failure  for  the  shear  events  in  the  district  are  (1)  right-lateral  strike-slip  movement 
along  steeply  dipping  northwest  striking  planes,  (2)  normal  faulting  along  moderate  dipping  northerly 
striking  planes,  and  (3)  thrust  faulting  along  randomly  striking,  shallow  dipping  planes.  P  and  T  axes 
derived  from  the  fault  plane  solutions  are  consistent  with  the  directions  of  principal  horizontal  stresses 
from  previous  in  situ  measurements  (Jung  et  al.,  1994). 


Figure  2.  Travel  time  curve  for  seismic  events  in  the  Coeur  d'Alene  district  based  on 
independently  known  hypocentral  locations. 
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Stoping  of  the  lead-zinc-silver  ore  in  the  district  takes  place  on  steeply  dipping  veins  at  an  average  depth 
of  1650  meters  below  surface.  The  country  rock  consists  primarily  of  stiff,  vitreous  quartzite  (Young's 
modulus  =  55  GPa,  uniaxial  compressive  strength  =  300  MPa).  There  are  also  smaller  pockets  of  sericitic 
quartzite  and  thin  bands  of  argillite  (both  of  which  are  considerably  softer  and  weaker  than  vitreous 
quartzite).  Scott  et  al.  (1993)  discussed  the  influence  that  these  geological  features  may  have  on 
rockbursting.  Based  on  overcoring  measurements,  observations  of  fracturing  around  raisebore  holes,  and 
natural  earthquake  data  the  maximum  principal  stress,  SI,  is  believed  to  lie  approximately  horizontal 
toward  N40W  (Whyatt  et  al.,  1993;  Sprenke  et  al.,  1991).  The  magnitude  of  SI  is  believed  to  be  on  the 
order  of  1.5  to  3.0  times  greater  than  the  vertical  stress.  The  magnitudes  and  orientations  of  the  other  two 
principal  stresses,  S2  and  S3,  are  less  certain. 

Three  distinct  types  of  seismic  monitoring  systems  are  used  to  monitor  rockburst  activity  in  the  mining 
district  (Williams  et  al,  1995).  Eiectrolab  MP250  microseismic  systems  were  installed  in  the  early  1970s 
to  monitor  rock  popping  around  stopes.  The  U.S.  Bureau  of  Mines  installed  in-mine  16-channel 
macroseismic  monitoring  systems  in  several  mines  since  1989.  This  mine-wide,  digital,  full-waveform 
system  provides  backup  source  locations  as  well  as  constraints  on  failure  mechanisms  and  peak-particle 
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Figure  3.  Longitudinal  section  view  of  the  Lucky  Friday  Mine  showing 
rockburst  hypocenters. 


ground  velocities.  The  University  of  Idaho  implemented  a  network  of  16  surface  seismographs  covering 
the  entire  Coeur  D'Alene  mining  district,  from  1992-1994.  All  the  surface  seismometers  were 
uncalibrated  vertical  component  velocity  gauges  ranging  in  natural  frequencies  from  1  Hz  to  4.5  Hz. 
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Figure  4.  Fault  plane  solutions  for  events  below  the  5400  sub-level  of  the  Lucky  Friday  Mine. 


This  network,  in  conjunction  with  a  surface  seismograph  in  the  mineyard,  provided  complete  focal-sphere 
coverage  of  seismic  events  at  the  Lucky  Friday  Mine,  and  allowed  the  computation  of  well-constrained 
fault  plane  solutions  for  most  events  above  magnitude  0.5. 

During  the  two-year  deployment  of  the  district-wide  seismograph  network,  the  mine  seismograph 
recorded  1750  seismic  events  with  hypocenter  locations  near  the  mine.  The  magnitudes  of  these  events 
ranged  from  -1  to  3.  The  events  were  assigned  locations  from  the  in-mine  rockburst  monitoring  systems. 
For  most  events,  the  source  location  is  accurate  to  within  30  meters  or  less.  Using  data  from  our  district¬ 
wide  seismograph  network,  fault  plane  solutions  were  attempted  for  each  of  the  1750  seismic  events 
recorded  on  the  surface  seismograph  at  the  Lucky  Friday  Mine  during  a  one-year  period.  Many  of  the 
events  proved  to  be  too  small  or  too  noisy  to  provide  impulsive  first  motions  on  the  far-field  stations  and 
some  events  resulted  in  ambiguous  or  poorly-constrained  solutions.  About  half  of  the  events  had 
implosional  source  components  (Stickney  and  Sprenke,  1993),  producing  dilatational  first  motions  in  all 
quadrants,  making  double-couple  fault  plane  solutions  impossible.  In  total,  well-constrained  fault  plane 
solutions  for  250  of  the  1750  events.  The  definition  of  "well-constrained"  was  that  there  were  no 
discrepant  stations  and  that  the  possible  solution  planes  fell  in  clusters  within  +/- 12  degrees  from  their 
mean  positions.  Fault  plane  solutions  were  calculated  by  the  program  FPFIT  (Lee  and  Valdes,  1985). 
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Figure  5.  Concentration  contours  for  the  P,  I,  and  T  axes  plotted  on  equal  area  stereonets. 

In  Figure  3,  the  locations  of  the  shear-type  seismic  events  are  projected  on  a  longitudinal  section  view  of 
the  bottom  portion  of  the  vein.  Fault  plane  solutions  corresponding  to  these  events  are  shown  in  plan 
view  (Fig.  4).  The  majority  of  the  events  occurred  near  the  three  bottom  longwall  stopes.  Two  types  of 
fault  plane  solutions  were  found  predominant:  (1)  randomly  oriented  normal  or  thrust  faulting,  and  (2) 
strike-slip  movement  along  vertical  planes  striking  N45W  or  N45E.  The  magnitude  3  event  located  near 
the  Offset  fault  shows  right-lateral  movement  on  a  plane  parallel  to  the  fault.  Also,  several  event  clusters 
adjacent  to  the  Offset  fault  all  show  the  same  strike-slip  movement,  located  along  a  straight  line  parallel 
to  die  fault. 

Because  virtually  all  of  the  250  shear  events  analyzed  occurred  in  the  vicinity  of  active  stopes,  the  effect 
of  near-  field  stress  concentrations  due  to  mine  geometry  was  critical  to  the  interpretation  of  each 
individual  event.  However,  a  general  analysis  of  the  fault  plane  solutions  in  conjunction  with  previous  in 
situ  measurements  of  stress  orientations  provided  some  interesting  observations  about  the  overall  stress 
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regime  in  the  mine.  The  respective  groupings  of  the  P,  I,  and  T  axes  for  all  250  well-constrained 
solutions  were  calculated  using  MicroNet  (Guth,  1987),  and  are  presented  in  lower  hemisphere,  equal 
area  stereonets  (Figure  5).  These  concentration  plots  indicate  that  the  preferred  seismic  P-axis  plunges 
near-  horizontally  with  a  bearing  about  N5W,  and  that  the  preferred  seismic  T-axis  is  near  horizontal, 
bearing  east-west.  In  terms  of  fault  planes,  this  indicates  the  preferred  failure  plane  is  vertical  and  strikes 
either  N50W  or  N40E.  In  terms  of  die  geology  of  the  mine,  die  most  reasonable  interpretation  is  that 
when  large  bursts  occur,  it  is  the  three  major,  northwest  trending,  strike-slip  faults  and  associated  parallel 
fractures  that  tend  to  be  re-activated  in  the  vicinity  of  active  stopes.  It  is  generally  easier  to  re-activate 
existing  rupture  planes  than  it  is  to  create  new  fractures.  The  old  fractures  tend  to  have  less  cohesion  and 
lower  angles  of  friction  than  the  surrounding  intact  rockmass.  If  the  slip  is  in  fact  on  pre-  existing 
northwesterly  trending  fractures,  the  fault  plane  solutions  show  that  the  movement  is  nearly  always  right- 
lateral  strike-slip  with  a  smaller  oblique  components  of  normal  or  reverse  faulting.  Furthermore,  because 
the  SI  direction  can  lie  anywhere  from  the  45°  position  of  the  seismic  P  axis  to  as  close  as  5°  to  the  fault 
plane,  the  preferred  P  and  T  axes  orientations  are  from  0°  to  40°  clockwise  from  the  actual  SI.  Previous 
in  situ  measurements  of  SI  in  and  around  the  mine  indicate  a  N40W  orientation  of  SI,  35°  from  the 
seismic  P  axis  and  10°  from  the  failure  plane,  in  the  range  one  might  expect  for  failure  on  pre-  existing 
faults. 

Composite  analysis  of  the  fault  plane  solution  data  and  previous  in  situ  stress  measurements  suggests  that 
the  overall  stress  regime  in  the  mine  is  as  follows:  horizontal  SI  oriented  N40W,  seismic  P  axis 
horizontal  N5W,  S2  vertical,  S3  oriented  horizontal  N50E,  and  seismic  T  axis  horizontal  N50E.  The 
overall  preferred  mode  of  failure  for  the  shear  events  at  the  Lucky  Friday  Mine  is  right-lateral  strike-slip 
moyement  along  steeply  dipping  N50W  striking  planes.  However,  significant  normal  faulting  along 
moderately  dipping,  northerly  striking  planes  and  thrust  faulting  along  randomly  striking,  shallow  dipping 
planes  also  occurs.  Except  for  a  few  normal  fault  events  on  the  5 100  level  near  the  North  Control  fault, 
most  of  the  fault  plane  solutions  do  not  correlate  well  with  the  bedding  plane  orientations.  Although  the 
direction  of  maximum  compression  sometimes  drifts  to  the  vertical  position  resulting  in  gravity  failures, 
the  east-to-northeast  extension  direction  is  particularly  persistent,  being  a  feature  of  most  of  the  fault 
plane  solutions. 

In  a  recent  study  (Whyatt  et  al.,  1996),  all  available  seismic  data  from  the  Lucky  Friday  Mine  for  the 
period  1989-1994  was  analyzed  to  provide  a  foundation  for  assessing  the  seismic  characteristics  of  large 
events  in  deep  mines.  Information  was  collected  and  analyzed  for  39  seismic  events  with  local 
magnitudes  above  2.5.  One  of  these  events  occurred,  on  average,  every  8  weeks  during  the  study  period. 
The  characteristic  first  motion  patterns,  damage  patterns,  and  relationships  to  mining  and  geologic 
structure  were  defined  for  each  event.  However,  other  than  the  first  motion  fault  plane  solutions,  no  focal 
mechanism  studies  or  moment  tensor  solutions  have  been  obtained  to  date  in  the  Coeur  d'Alene  district. 
Furthermore,  no  high  quality  broadband  instruments  have  ever  been  deployed  in  or  near  the  district. 
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1.4  Objectives. 

The  objectives  of  this  phase  of  the  contract  were 

1.  To  deploy  a  local  network  of  seismic  sensors  in  a  high  rockburst  region  near  existing  mines  in 
Northern  Idaho. 

2.  To  collect  seismic  data  from  this  sensor  network,  from  in-mine  monitoring  stations,  and  from  regional 
and  national  networks  at  distances  from  tens  to  hundreds  of  kilometers. 


Figure  6.  Vertical  Waveforms  from  Rockbursts  Associated  with  Blasting  at  the  Lucky  Friday  Mine. 


Locations  and  focal  mechanisms  determined  from  in-mine  and  regional  seismic  arrays  indicate  a 
significant  implosional  component  (in  addition  to  the  usual  shear  component  of  natural  earthquakes)  for 
rockbursts  and  mine  tremors  in  the  Coeur  d’Alene  mining  district  (Figure  1)  of  northern  Idaho  (Sprenke  et 
al.,  1991;  Stickney  and  Sprenke,  1993;  Jung  et  al.,  1995).  A  small,  non-double-couple  component 
seismic  event  at  shallow  depth  (less  than  a  few  kilometers)  may  be  difficult  to  discriminate  from  an 
explosion  at  regional  distances.  The  initial  P-wave  polarity  may  not  be  determinable,  the  location  and 
waveform  characteristics  will  indicate  a  source  in  the  shallow  crust,  and  the  relative  amplitudes  of 
compressional  and  shear  waves  may  be  similar  to  that  of  an  explosion.  There  have  been  an  average  of 
approximately  two  seismic  events  per  year  since  1990  (Table  1)  from  the  Lucky  Friday  mine  that  have 
magnitudes  larger  than  3  and  these  events  could  be  detected  by  CTBT  monitoring  systems. 
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Table  1.  Event  catalog  from  the  Composite  National  Seismic  System 


Date 

Time 

Lat 

Long 

Depth 

Mag 

3/7/1990 

21:54.4 

47.453 

-115.974 

1.7 

3 

6/13/1991 

03:36.7 

47.852 

-115.713 

7.5 

3 

9/19/1991 

20:27.8 

47.458 

-115.803 

1 

3 

12/11/1991 

37:03.4 

47.48 

-115.801 

1 

3 

8/4/1992 

26:28.9 

47.468 

-115.785 

2 

3.2 

8/12/1992 

30:23.4 

47.469 

-115.78 

2.3 

3.2 

11/4/1992 

21:42.7 

47.467 

-115.775 

2.2 

3 

3/18/1993 

00:53.1 

47.472 

-116.015 

9.6 

3.1 

wnaEmm 

in 

o 

o 

00 

47.477 

-116.108 

1 

3.2 

10/18/1993 

49:10.1 

47.471 

-115.783 

1 

3.3 

3/8/1994 

07:11.5 

47.186 

-115.932 

5 

3.4 

8/17/1994 

56:49.0 

47.472 

-115.778 

1.6 

4.1 

4/25/1996 

56:29.1 

47.443 

-115.816 

1 

3 

5/30/1996 

00:37.8 

47.44 

-115.82 

1 

3.6 

11/2/1996 

28:26.3 

47.798 

-115.642 

5 

3.6 

3/14/1997 

18:32.0 

47.47 

-115.78 

2 

3.6 

This  project  obtained  the  first  high-quality  seismic  data  to  determine  the  source  characteristics  and 
moment  tensors  of  seismic  events  in  this  region.  This  is  the  first  time  a  broad-band,  high  dynamic-range, 
three-component  instrument  has  been  operated  in  this  mining  district. 

Table  2.  Event  Catalog  from  the  In-Mine  Monitoring  System 


DATE 

TIME 

X 

Y 

Z 

AMPL  (mm) 

STOPE 

ML 

7/30/1998 

1:20:28 

232 

-164 

-2750 

7 

575-05 

7/30/1998 

16:24:12 

287 

-1000 

-2564 

4 

550-01 

7/30/1998 

23:40:13 

-397 

-100 

-2750 

3 

575-05 

7/31/1998 

1:33:09 

-356 

-307 

-2750 

4 

575-05 

7/31/1998 

5:25:36 

-396 

-1000 

2750 

10 

550-01 

7/31/1998 

23:35:54 

-988 

-460 

-2750 

4 

550-01 

8/1/1998 

5:53:22 

-258 

-71 

-2715 

3 

575-05 

8/2/1998 

1:52:25 

-220 

-90 

-1463 

6 

49-87 

8/3/1998 

9:12:44 

-218 

-736 

-2760 

2 

550-01 

8/4/1998 

20:44:42 

-92 

-763 

-2750 

3 

550-01 

8/4/1998 

21:29:38 

1793 

108 

-1380 

5 

46-110 

8/4/1998 

22:18:09 

-663 

-320 

-2750 

2 

550-01 

8/5/1998 

0:39:29 

-177 

723 

-2750 

14 

593 -ramp 

8/5/1998 

23:31:47 

-169 

587 

-2750 

4 

593 -ramp 

8/6/1998 

7:25:31 

-54 

-512 

-2585 

74 

575-05 

mw+m 

8/6/1998 

18:19:23 

156 

-1000 

-2750 

5 

550-01 

B 

8/8/1998 

1:54:56 

-1048 

2352 

-1726 

76 

paymaster 

8/10/1998 

23:37:42 

-255 

496 

-2750 

3 

5 93 -ramp 

8/10/1998 

23:38:54 

220 

-1 66 

-2750 

2 

584-06 
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Table  2.  Event  catalog  from  the  in-mlne  monitoring  system  (eont.) 


DATS 


a/11/1998 


8/12/1998 


8/13/1998 


8/13/1998 


8/13/1998 


8/14/199 


8/16/1998 


8/18/199 


8/18/1998 


/18/199 


8/19/1998 


/20/199 


8/20/1998 


/20/199 


8/2Q/1998 


/21/199 


8/21/1998 


/21/199 


8/21/1998 


8/22/199 


8/22/1998 


/ 24/1998 


8/24/1998 


8/25/199 


8/25/1998 


8/26/199 


8/26/1998 


8/26/1998 


8/27/1998 


8/27/1998 


8/27/1998 


8/28/1998 


8/28/1998 


8/28/1998 


8/28/1998 


8/28/1998 


8/28/1998 


8/28/1998 


8/29/1998 


11:27:33 


13:32:46 


3:59:22 


8:54:49 


14:35:49 


2:54:31 


1:05:42 


11:25:28 


17:18:34 


23:39:14 


22:51:55 


13:42:43 


13:42:57 


14:09:26 


20:18 


31:35:00 


31:38:00 


34:23:00 


36:47:00 


38:26:00 


33:59:00 


11:12:13 


16:04:47 


13:35:19 


23:36:01 


0:18:48 


13:55:57 


14:12:32 


0:08:09 


0:57:41 


14:47:24 


0:24:03 


7:54:22 


13:35:40 


13:49:50 


13:53:29 


19:13:15 


23:09:27 


23:24:16 


23:27:41 


7:26:56 


9:07:36 


13:56:06 


18:57:40 


11:53:25 


-6 


-33 


-729 


1146 


-32 


-1445 


-255 


183 


-281 


-153 


-161 


-248 


-106 


93 


124 


-137 


152 


139 


-335 


-323 


1985 


511 


-342 


-200 


165 


-235 


-193 


-154 


-320 


-130 


-42 


-160 


203 


19 


83 


-82 


53 


357 


98 


50 


-216 


1514 


-600 


14 


Y 


-205 


197 


-415 


70 


-312 


3150 


-980 


-162 


-70 


-47 


-438 


-107 


-146 


87 


-261 


-286 


-228 


-66 


-890 


-306 


-200 


1605 


299 


-168 


-195 


-140 


-247 


-168 


-311 


-220 


-195 


-36 


-9 


-307 


46 


-347 


-73 


81 


-65 


-237 


-432 


-37 


29 


55 


-2750 


-2750 


-2750 


-1989 


-2750 


-1638 


-2423 


-2433 


-2750 


-2323 


-2750 


-2750 


-2750 


-2750 


-2378 


-2419 


-2749 


-2750 


-2750 


-2750 


-2750 


-1322 


-2584 


-2750 


-2752 


-2750 


-2731 


-2716 


-2750 


-2750 


-2738 


-2750 


-2081 


-2750 


-2547 


-2348 


-2750 


-2075 


-2750 


-2750 


-2750 


-2750 


-2750 


-2750 


-2309 


575-05 


584-06 


5900 


575-05 


1.5  Seismic  Data  from  Sensor  Network  and  from  In-mine  Stations. 

We  operated  five  strong-motion  accelerometers  and  one  broadband  seismometer  within  about  5  km  of  the 
Lucky  Friday  mine  during  the  summer  of  1998.  The  strong  motion  accelerometers  recorded  2-10  of  the 
larger  rockbursts  per  week  of  operation;  the  broadband  instrument  recorded  about  20-40  per  week.  We 
have  collected  about  24  weeks  of  data  from  the  broadband  seismometer  and  about  14  weeks  of  data  from 

*  the  strong  motion  accelerometers  (Table  2).  Many  of  the  events  occured  at  or  near  the  time  explosions 
were  used  to  create  new  openings  in  the  mine  (Table  3),  but  most  waveforms  indicate  implosional  first 
motions.  The  explosions  were  usually  too  small  to  be  observed,  and  the  implosional  events  often  occured 

*  several  minutes  after  the  blasting  time.  Figure  6  shows  examples  of  the  waveforms  obtained  on  the 
broadband  system  for  events  occurring  near  the  blast  times.  Figure  7  shows  the  displacement  spectra 
obtained  from  26  events  that  all  had  similar  waveforms  and  dilational  first  motions  at  this  site.  The  use  of 
broad-band,  high  dynamic  range  accelerometers  produced  very  high  quality,  on-scale  recordings  of 
rockbursts  at  distances  as  short  as  2  km  and  magnitudes  up  to  2. 

The  Lucky  Friday  survey  system  is  (X,Y,Z)  where  X  =  Easting  (ft),  Y  =  Northing  (ft),  and  Z  = 
Elevation  (ft).  Conversion  from  mine  coordinates  to  geodetic  survey  coordinates  is  described  in  Table  4. 
Some  mine  maps  will  show  coordinates  with  20000  added  to  X  and  Y,  so  the  origin  used  by  the  MP250  is 
technically  at  (20000,20000,0)  on  these  maps.The  X  and  Y  used  in  the  conversions  below  should  NOT 
have  the  20000  added.  Note  that  most  bursts  will  be  at  negative  elevations.  The  elevation  of  the  surface 
seismograph  at  the  mine  is  3360  feet  above  sea  level. 

Table  3:  Blast  times  for  mines  in  the  Coeur  d'Alene  District 


Gold  Hunter  Mine 

115  am  PDT 

0815  UTC 

130  am  PDT 

0830  UTC 

Gold  Hunter  Mine 

1100  am  PDT 

1800  UTC 

Gold  Hunter  Mine 

2215  UTC 

E5B5E5M 

2230  UTC 

Gold  Hunter  Mine 

Easaasai 

0400  UTC 

Table  4:  Mine  coordinates 


LATITUDE 

0.0001638*Y  +  47  degrees  28.18  minutes 

LONGITUDE 

115  degrees  46.90  minutes  -  0.0002427*X  +  0.000004*Y 

ELEVATION  (meters) 

Z  *  0.3048 

A 


V* 


11 


frequency.  Hi 


Figure  8.  Waveforms  recorded  at  the  Yellowknife  Array,  1500  km  north  of  the  Lucky  Friday  Mine. 
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1.6  Data  from  Sources  at  Regional  Distances. 

We  have  also  obtained  some  regional  waveforms  from  North  American  stations  that  are  elements  or 
possibly  auxiliary  elements  of  the  CTBT  monitoring  system.  Figure  8  shows  the  waveforms  recorded  at 
the  Yellowknife  array  from  an  event  at  die  Lucky  Friday  mine  in  1996.  Data  from  other  North  American 
broadband  stations,  and  regional  arrays  operated  in  eastern  Washington,  northern  Idaho,  and  western 
Montana,  were  assembled  to  evaluate  wave  propagation  characteristics  and  the  development  of  regional 
phases  from  seismic  events  in  the  study  region.  A  catalog  of  all  known  recorded  rockbursts  in  the 
northern  Idaho  region  is  given  in  Table  5  below. 

The  Montana  Bureau  of  Mines  and  Geology  (MBMG)  records  a  network  of  seismograph  stations 
operated  in  southwest  Montana  and  collects  trace  data  from  two  smaller  networks  operated  in  west-central 
and  northwest  Montana.  Together  these  three  groups  of  seismograph  stations  comprise  the  Montana 
regional  seismograph  network.  Local  and  regional  earthquakes  were  visually  identified  on  seismograms 
and  P  and  S  arrival  time  were  picked  if  three  or  more  stations  record  first-arriving  P-waves.  Mining  and 
construction  blasts  were  discriminated  on  the  basis  of  seismogram  appearance,  time  of  day,  and  source 
proximity  to  known  mining  or  construction  areas.  Hypocenter  locations  were  determined  for  all  seismic 
events  that  could  not  be  confidently  identified  as  blasts  and  are  included  in  the  MBMG  earthquake 
catalog.  Phase  data  were  collected  from  seismograph  stations  operated  in  surrounding  areas  by  other 
institutions.  Surrounding  seismograph  stations  include  those  operated  by  the  US  National  Network 
Station  at  Newport,  WA,  the  University  of  Washington,  the  University  of  Idaho,  the  Army  Corps  of 
Engineers,  Boise  State  University,  the  Idaho  National  Engineering  and  Environmental  Laboratories, 
surface  mine  seismographs  at  three  mines  in  the  Coeur  d'Alene  mining  district  operated  by  various  mining 
companies  and  the  U.S.  Bureau  of  Mines,  and  various  stations  in  southern  British  Columbia  and  southern 
Alberta  operated  by  the  Canadian  Geological  Survey.  Of  course  not  all  of  these  stations  operated  during 
the  same  time  periods.  Seismograph  coverage  in  the  Northern  Rocky  Mountain  region  has  generally 
improved  over  the  past  18  years  however,  the  best  seismic  coverage  of  the  Coeur  d'Alene  mining  district 
existed  while  the  North  Idaho  seismic  network  was  operational  in  the  early  1990s. 

P  and  S  arrival  times  were  initially  read  from  paper  seismograms  recorded  at  a  rate  of  60  mm/second  for 
an  overall  timing  accuracy  of  0.1  seconds.  After  August  of  1989,  Montana  network  stations  were  digitized 
at  100  samples/second  and  recorded  on  a  PC-based  data  acquisition  in  Butte.  XRTP  (Tottingham,  1994), 
file  data  acquisition  software  operated  in  a  triggered  mode.  Events  that  triggered  were  picked  using  PC- 
SUDS  software  (Banfill,  1996)  for  an  overall  timing  accuracy  of  0.02  seconds  for  impulsive  arrivals. 
When  the  Flathead  network  became  operational  in  1996,  XRTP  was  used  to  record  seismic  events.  XRTP 
was  also  installed  in  Missoula  at  about  this  time  to  record  west-central  Montana  seismic  data.  In  1998, 
both  recording  sites  were  upgraded  to  ViSeis  data  acquisition  software.  Phase  data  from  all  available 
seismograph  stations  were  used  with  the  Montana  crustal  velocity  model  (Stickney,  1997)  and  the 
hypocenter  location  programs  HYP071  (Lee  and  Lahr,  1975)  and  HYP071PC  (Lee  and  Valdes,  1985)  to 
determine  hypocenter  locations.  Magnitudes  were  determined  from  coda  length  measurements  using  the 
method  described  by  Stickney  (1995).  When  visible  signals  were  recorded  on  the  Butte  Wood-Anderson 
equivalent  seismograph,  local  (Richter)  magnitudes  were  also  determined.  Magnitudes  determined  from 
both  methods  for  individual  events  generally  agree  within  one  quarter  of  a  magnitude  unit. 

A  search  of  the  MBMG  earthquake  catalog  of  an  area  10  km  in  radius  around  the  Lucky  Friday  Mine 
(47.471°N,  1 15.783°W)  yielded  86  seismic  events  (Table  5).  All  are  believed  to  be  seismic  events  directly 
associated  with  the  Lucky  Friday  Mine.  Mine  operators  confirmed  that  many  of  these  events  were 
rockbursts.  Of  these  86  events,  two  resulted  in  causalities  to  underground  miners,  and  19  were  reported  as 
being  felt  by  people  at  the  surface,  although  many  more  were  undoubtedly  felt  but  not  so  reported.  The 
cataloged  Lucky  Friday  Mine  seismic  events  since  1982  range  in  magnitude  from  1.5  to  4.1  although  the 
detection  threshold  for  the  Montana  regional  seismograph  network  alone  is  about  magnitude  2.3.  A  total 
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of  18  events  have  magnitudes  of  3.0  or  greater,  6  are  3.5  or  greater  and  1  exceeds  magnitude  4.0.  From 
1996  through  1999,  on  average,  8.8  Lucky  Friday  events  are  detected  and  located  annually  with  regional 
network  data. 


Table  5.  Regional  Rockburst  Catalog  (1982-1999) 


Date 

hr-min 

sec 

Latitude 

Depth (km) 

Magnitude 

83 

10 

15 

1115 

10.3 

47,485 

115.79 

0.2 

2,5 

84 

1 

19 

725 

47.1 

47.401 

115.839 

0.4 

3,1 

84 

1 

25 

1936 

11 

47.488 

115.734 

11 

2.9 

85 

1 

20 

118 

24.6 

47.526 

115.767 

1 

2.8 

85 

2 

8 

1422 

19 

47.469 

115.81 

4.4 

3.2 

85 

4 

13 

510 

31 

47.439 

115.817 

0.1 

1.7 

85 

9 

12 

2134 

40.5 

47.484 

115,775 

9.5 

1.5 

85 

10 

31 

1131 

50.9 

47.514 

115.751 

4.7 

3 

85 

12 

4 

2340 

51.9 

47.494 

115.776 

6.3 

3.2 

85 

12 

10 

1858 

43.5 

47.413 

115.867 

0.2 

2.7 

86 

'  1 

9 

634 

53.8 

47.477 

115.771 

6.4 

2.9 

86 

1 

28 

1324 

0.3 

47.486 

115.722 

7.5 

2.8 

86 

1 

30 

2239 

26.4 

47.506 

115.722 

5.7 

2.8 

86 

3 

5 

1939 

42.7 

47.483 

115.749 

0.6 

2.1 

86 

3 

12 

1632 

56.7 

47.46 

115.802 

0.1 

2.6 

86 

3 

25 

2207 

38.8 

47.472 

115.746 

1.3 

2.3 

86 

5 

13 

506 

6 

47.453 

115.9 

1.8 

2.5 

87 

11 

19 

637 

29.6 

47.456 

115.838 

2.6 

2.6 

87 

11 

28 

321 

35.6 

47.499 

115.772 

14.8 

2.4 

88 

2 

3 

1506 

52.5 

47.427 

115.747 

1.2 

2.4 

88 

10 

7 

540 

27.4 

47.428 

115.887 

10.6 

2.5 

88 

10 

18 

2115 

1.7 

47.463 

115.826 

8.3 

3.7 

89 

4 

4 

1042 

59,3 

47.462 

115.846 

7.1 

3.1 

89 

4 

13 

1323 

32.8 

47.464 

115.853 

9 

2.5 

89 

7 

6 

707 

56 

47.462 

115.839 

8.9 

3.7 

89 

10 

16 

1101 

7 

47.472 

115.863 

9,3 

2.6 

89 

11 

3 

631 

1.1 

47.459 

115.814 

4.2 

2.8 

90 

3 

31 

703 

5.7 

47.425 

115.84 

3.5 

2.4 

90 

4 

11 

422 

42.9 

47.468 

115.77 

0.9 

2.6 

90 

6 

12 

2104 

56.2 

47.455 

115.804 

2.9 

2.5 

90 

8 

3 

1639 

13.8 

47,471 

115.783 

1,5 

2.8 

91 

2 

23 

652 

13.1 

47.434 

115.853 

13.9 

3 

91 

3 

1 

409 

25,6 

47.42 

115.694 

14.9 

2.8 

91 

5 

9 

1918 

7 

47.472 

115.842 

8.7 

2.8 

91 

5 

17 

734 

16.4 

47,45 

115.806 

3 

2.7 

91 

8 

17 

832 

51.3 

47.452 

115.892 

9,7 

2.9 

91 

9 

19 

1620 

28.8 

47.455 

115.826 

3.4 

3.5 

91 

12 

11 

2137 

0 

47.5 

115.8 

2 

3.2 

92 

7 

30 

738 

20.8 

47.468 

115.785 

2.2 

2.8 

92 

8 

4 

926 

28.9 

47.468 

115.785 

2 

3.2 

92 

8 

12 

830 

23.4 

47.469 

115.78 

2.3 

3.2 

92 

11 

4 

521 

42.7 

47.467 

115.775 

2.2 

3 

93 

10 

18 

649 

10.1 

47.471 

115.783 

1 

3.3 

94 

8 

17 

656 

50.3 

47.472 

115,778 

1.6 

4.2 

94 

9 

28 

342 

7 

47.471 

115.783 

2 

2.8 

94 

12 

4 

455 

59.8 

47.471 

115.77 

2 

2.8 
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Table  5.  Regional  Rockburst  Catalog  (1982-1999)  (Continued) 


Date 


95  1  13 


95  4  26 


95  5  27 


95  6  23 


95  9  19 


96  2  9 


96  4  25 


96  5  30 


96  9  26 


96  9  26 


96  12  20 


97  2  14 


97  3  14 


97  . 3  27 


97  4  26 


97  5  25 


97  6  28 


97  8  22 


97  9  5 


97  9  24 


97  11  8 


97  11  12 


98  1  30 


98  5 


98  5  14 


98  8  29 


98  12  9 


98  12  18 


99  1  25 


99  6  12 


99  10  23 


99  10  28 


99  11 


99  12  17 


hr~min 


836 


1414 


1117 


438 


2040 


201 


2056 


1000 


2137 


2259 


240 


838 


837 


1718 


2109 


558 


1534 


302 


2128 


2146 


2250 


1333 


2314 


840 


2336 


1811 


937 


1200 


2237 


834 


609 


2207 


2316 


2214 


1843 


2042 


2216 


707 


847 


1025 


Latitude 

Longitude 

47.471 

115.783 

47.471 

115.783 

47.471 

115.783 

47.44 

115.855 

47.457 


47.445 


47.471 


47.471 


47.471 


47.468 


47.471 


47.471 


47.471 


47.46 


47.542 


47.471 


47.471 


47.471 


47.471 


47.471 


47.415 


47.48 


47.471 


47.471 


47.471 


47.471 


47.471 


47.471 


47.471 


47.454 


47.471 


47.471 


47.471 


47.471 


47.471 


47.475 


47.469 


47.489 


47.436 


115.783 


115.84 


115.807 


115.783 


115.783 


115.783 


115.845 


115.783 


115.783 


115.783 


115.864 


115.726 


115.783 


115.783 


115.783 


115.783 


115.783 


115.818 


115.763 


115.783 


115.783 


115.783 


115.783 


115.783 


115.783 


115.783 


115.833 


115.783 


115.783 


115.783 


115.783 


115.783 


115.805 


115.806 


115.777 


115.827 


Depth (km) 


2 


2 


2 


2 


2 


0.2 


.3 


2 


2 


2 


2 


2 


2 


2 


2 


4.5 


2 


2 


2 


2 


2 


.5 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


Magnitude 


2.8 


SECTION  2 


SEISMIC  SOURCE  PARAMETERS  OF  ROCK  BURSTS 
COEUR  D’ALENE  MINING  DISTRICT,  IDAHO 


2.1  Introduction. 

Characterization  of  mining-related  seismic  events  is  critical  to  the  implementation  of  the  Comprehensive 
Nuclear  Test-Ban  Treaty  monitoring  system.  Mine  seismicity  from  uncontrolled  sources  are  particularly 
troublesome,  with  magnitudes  up  to  4  or  5,  producing  signal  strengths  comparable  to  1-10  kiloton 
contained  nuclear  explosions.  In  this  study,  ground  motions  from  four  surface  stations  and  one 
underground  station  (Figure  9)  were  inverted  to  determine  the  complete  moment  tensor  of  a  relatively 
large  rock  burst  in  the  Coeur  d’Alene  Mining  District  in  northern  Idaho.  Many  rock  bursts,  induced  by  the 
mining  process,  show  shear  implosional  source  mechanisms,  a  feature  that  should  be  useful  to 
discriminate  mining-related  events  from  underground  explosions. 


Figure  9.  Locations  of  stations  and  of  larger  events  recorded. 

2.2  Objective. 

The  objective  of  this  portion  of  the  study  was  to  apply  modem  methods  and  techniques  to  determine  the 
seismic  source  parameters  of  rockbursts.  Mining  activities  world-wide  produce  seismic  events,  some  of 
which  may  be  visible  to  monitoring  for  the  Comprehensive  Test  Ban  Treaty  (CTBT).  Mine  seismicity 
from  uncontrolled  sources  are  particularly  troublesome,  with  magnitudes  as  high  as  4  to  5,  producing 
signals  equivalent  to  1-10  kiloton  contained  nuclear  explosions.  Rock  bursts  in  deep  underground  mines 
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are  in  this  latter  category.  These  events,  however,  may  have  characteristics  unlike  either  explosions  or 
earthquakes  that  might  be  used  to  “fingerprint”  specific  mines  for  discrimination  purposes.  In  this  light, 
we  are  investigating  seismic  signals  from  one  of  the  most  rock  burst-prone  hard  rock  mining  areas  in  the 
world:  the  Coeur  d’Alene  Mining  District  of  Northern  Idaho. 
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Figure  10,  Displacements  recorded  at  MOR  for  the  largest  rockburst  recorded. 


In  a  review  of  mining-induced  seismicity,  Gibowicz  and  Kijko  (1994)  note  that  such  events  generally  fall 
into  two  types — those  directly  associated  with  mining  operations  and  those  associated  with  movement  on 
major  geological  discontinuities.  Rock  bursts  in  the  Coeur  d’Alene  District  are  not  so  easily  categorized. 
While  the  events  clearly  are  associated  with  mining  rate  (Sprenke  and  others,  1991)  and  stope  geometry 
(Jung  et  al.,  1995),  they  have  long  been  known  to  be  the  result  of  failure  on  faults  and  bedding  planes 
(Scott,  1990)  associated  with  the  Osbum  Fault,  part  of  the  Lewis  and  Clark  Line  megashear  that  strikes 
north-northwest  through  the  district.  An  underhand-Iongwall  system  of  mining,  implemented  in  the  late 
1980s  at  the  Lucky  Friday  Mine,  has  greatly  reduced  rock  burst  problems  at  that  mine,  illustrating  the 
influence  of  mining  geometry  on  rock  bursts  .  A  district-wide  network  operated  in  the  early  1990s  (Jung 
et  al.,  1995  )  showed  that  those  shear  events  that  could  be  fit  by  a  double  couple  solution  were  consistent 
with  a  regional  tectonic  stress  with  a  northwesterly  oriented  seismic  P-axis  and  and  northeasterly  seismic 
T-axis,  an  orientation  particularly  deleterious  for  slip  on  the  north-northwesterly  oriented  vertical  faults 
and  steeply-dipping  bedding  parallel  to  the  Osbum  Fault.  Apparently,  mining  operations  in  the  Coeur 
d’Alene  district  are  releasing  pent-up  strain,  strongly  affected  by  local  geology  and  tectonics.  Not  unlike 
deep  mines  world-wide,  the  seismicity  in  this  district  are  affected  by  depth  of  mining.,  production  rate, 
mining  geometry,  geologic  structure,  and  geologic  discontinuities. 

Although  blasting  often  triggers  rock  bursts  instantaneously,  many  rock  bursts,  particularly  large 
damaging  events  in  deep  underground  mines,  occur  spontaneously  and  unexpectedly — this  despite  the 
fact  that  they  are  caused  by  stress  changes  associated  with  file  mining  process.  This  interaction  of  mine 
geometry  and  geological  structure  suggests  that  rock  bursts  in  deep  underground  mines  may  often  be  of  a 
shear-implosional  nature,  involving  both  co-seismic  opening  closure  and  shear  failure  in  the  surrounding 
strata  (e.g.  McGarr,  1992).  Seismograph  networks  (Stickney  and  Sprenke,  1992;  Jung  et  al.,  1995) 
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demonstrated  that  well-constrained  all-dilatational  P-wave  first  motions  for  about  half  of  the  events  in  the 
Lucky  Friday  Mine  indicated  a  significant  implosional  component. 

For  this  project,  we  have  investigated  the  seismic  signals  from  a  damaging  rock  burst  in  the  Lucky  Friday 
Mine.  The  local  magnitude  of  the  event  was  3.1  as  measured  at  Butte,  Montana.  It  was  recorded  on  the 
Montana  and  Washington  region  networks.  Predominantly  dilatational  first  motions  for  this  event  from 
our  five  local  stations  and  several  regional  stations  suggested  that  this  event  had  an  implosional 
component.  The  following  source  parameters  were  estimated  from  the  events  far  field  displacement 
spectrum  as  measured  on  the  broadband  instrument  at  station  MOR  (Figures  9  and  10). 


23  Seismic  Moment. 

The  strength  of  a  seismic  event  is  commonly  defined  in  terms  of  a  double-couple  point  source  model. 
More  complex  source  characteristics  will  be  considered  later  in  this  paper  when  the  seismic  moment 
tensor  is  computed  for  the  events.  The  seismic  moment  Mo  is  the  usual  parameter  used  to  express  the  size 
of  a  seismic  event.  Moment  is  defined  by  Mo  =  p  u  A  where  p.  is  the  shear  modulus,  u  is  the  average  slip 
of  the  fault,  and  A  is  the  area  of  the  fault.  For  large  events,  M0  can  be  calculated  directly  from  field 
estimates  of  fault  displacement  and  fault  area,  but  usually  Mo  has  to  be  estimated  using  seismic  records  in 
the  frequency  domain  using  the  following  formula: 


Mo  —  4  7i  c3  R  Q0  Fc  1  Rc  *  Sc 1  (2.1) 


where 

c  =  velocity  of  P  waves  (a)  or  S  waves  ((3)  at  the  source, 

R  =  hypocentral  distance,  Q0  is  the  low  frequency  limit, 

Fc  =  the  radiation  factor, 

Rc  =  the  free  surface  effect, 

Sc  =  the  site  correction.  For  the  broadband  station  MOR  (Figure  9)  used  in  this  study, 

R  =  2660  m, 
a  —  5740  m/s, 

P  =  3470  m/s, 

Fca  =  0.52, 

FcP  =  0.63, 

RCa  =  1.05,  and 
RCp  =  2. 

The  low  frequency  limit  fio  can  be  measured  directly  on  the  far  field  displacement  spectrum  or  can  be 
found  by  integrating  the  displacement  pulse  in  the  time  domain.  For  our  largest  event,  Qo  for  P  was  found 
t°  be  1-1  x  10  m-s  and  for  S  was  1 .4  x  1C6  m-s.  The  corresponding  seismic  moments  were  calculated  as 
3.4  xlO  N-m  and  4.2  x  1012  N-m  respectively.  For  comparison,  as  described  later,  the  moment  tensor 
estimate  of  Mo  is  about  5  x  10u  N-m. 
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2.4  Moment  Magnitude. 


Moment  magnitude  Mw  is  commonly  used  as  a  measure  of  seismic  event  size.  The  definition  of  this 
magnitude  scale  is  as  follows: 


Mw=  2/3  log  Mo -6.07 


(2.2) 


where  Mo  is  expressed  in  N-m  (Kanamori  and  Hanks,  1979).  For  the  largest  rockburst  in  our  study,  Mw 
was  calculated  from  our  values  of  M0  as  2.4.  For  comparison,  the  magnitude  was  measured  on  the  Richter 
scale  as  ML  2.9  (Butte). 


2 £  P-Wave  and  S-Wave  Energy. 

The  radiated  elastic  energy  Ec  by  P  or  S-waves  is  another  measure  of  the  size  of  a  seismic  event.  The 
radiated  energy  can  be  approximated  as  follows: 

Ec = 4  n  p  c  <FC>2  R2  Fc'2  Rc*2  Jc  (2.3) 


where 

c  =  a  or  p,  and 

<FC>2  =  the  mean  square  radiation  coefficients,  4/15  for  a  or  2/5  for  p. 
The  other  parameters  are  as  defined  above,  and: 


Jc  =  0.25  Qo2  (2nf0)3  (2.4) 

where  f0  -  the  comer  frequency  of  the  displacement  spectrum. 

For  die  largest  event  in  our  data  set,  we  found  fo  —  20  Hz  for  the  P  wave  and  14.8  Hz  for  the  S  wave. 
Using  other  parameters  as  given  above,  the  P  and  S  radiated  energy  were  found  to  be  1500  N-m  and 
39,000  N-m,  respectively.  The  ratio  of  Es  to  EP  for  the  event  was  about  26,  quite  a  high  ratio  considering 
the  implosional  nature  of  the  source.  Rockbursts  are  said  to  be  commonly  depleted  in  S  waves  with  Es  to 
Ep  ratios  as  low  as  1.5  (Gibowicz  et  al.,  1990). 


2.6  Radiation  Efficiency. 

Rockbursts  in  mines  are  thought  to  have  low  seismic  efficiency.  Only  a  few  per  cent,  at  most,  of  the  • 

rupture  energy  seems  to  be  released  by  seismic  waves  (McGarr,  1976;  McGarr  et  al.,  1979).  Most  of  the 
energy  is  lost  as  heat.  The  radiation  efficiency  eC  is  a  measure  comparing  radiated  seismic  energy  for  P  or 
S  waves  to  dynamic  stress  drop  (srms )  at  the  source: 

Sc  =  pEcMo/arms  (2S) 
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where  p.  =  shear  modulus. 

The  dynamic  stress  drop  is  calculated  from  aRMS,  the  root  mean  square  acceleration  averaged  over  the 
duration  of  the  S  wave  arrival  as  follows: 


<Jrms  -  3.2  p  R  aRMs  (fs  /  f  max)1/2 


(2.6) 


where 

fs  =  S  wave  comer  frequency  on  the  displacement  sprectrum,  and 
fjnAX  =  highest  frequency  of  acceleration  reliably  recorded. 

For  the  largest  rockburst  recorded  in  our  study,  we  found  aRMS  =0.11  m/s2,  fs  =  14.8  Hz,  fMAX=  80  Hz 
and  Orms  =  1.08  x  106Pa.  The  radiation  efficiency  for  P  and  S  waves  was  found  to  be  1.33  x  10'5  and 
2.82  x  10"4,  respectively. 


2.7  Source  Dimensions. 

Estimates  of  source  dimension  are  commonly  made  assuming  the  Brune  (1970)  model,  for  which  the  area 
of  a  circular  fault  is  inversely  proportional  to  S  wave  comer  frequency: 


2n  r0  =  2.34  $/f0 


(2.7) 


where  r0  is  the  radius  of  the  fault.  This  relation  assumes  an  instantaneous  stress  release. 

For  the  largest  rockburst  in  our  data  set,  the  source  size  worked  out  to  be  87  m  using  the  simple  Brune 
model.  It  must  be  noted,  however,  as  McGarr  (1984, 1991)  points  out,  that,  although  displacement  may 
occur  over  the  entire  source  dimension,  the  actual  damge  in  mines  would  occur  over  a  much  smaller  area 
where  asperities  along  the  fault  exist. 


2.8  Pulse  Duration  and  Rupture  Velocity. 

An  estimate  of  the  upper  limit  of  rupture  velocity  can  be  made  from  x,  the  pulse  duration  of  ground 
velocity. 


V=  12r0/  (16  x  +  13  £)  (2.8) 

where  4  =  r0  sin  0  /  p. 

The  angle  0  is  the  take  off  angle  measured  from  the  normal  to  the  fault  plane.  If  0  =  Oo,  then  the  upper 
limit  of  rupture  velocity  can  be  estimated.  For  our  largest  rockburst,  the  takeoff  angle  0  is  close  to  0°.  The 
rupture  velocity  comes  out  at  1770  m/s  or  about  0.51  p,  somewhat  slower  than  natural  earthquakes  which 
have  rapture  velocities  from  0.6  p  to  0.9  p. 
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2.9  Stress  Drop. 

The  difference  between  the  initial  and  final  stress  levels  over  a  fault  plane  is  called  the  static  stress  drop. 
However,  seismic  data  can  really  only  measure  the  dynamic  stress  drop,  the  difference  between  the  initial 
stress  level  and  the  kinetic  friction  level  on  the  fault. 

The  most  common  estimate  of  stress  release  during  a  seismic  event  is  the  Brune  Stress  Drop  A  a 
calculated  from: 


Act  =  7/16  M0/r03  (2.9) 

The  Brune  Stress  Drop  is  an  approximation  of  the  static  stress  drop  and  it  represents  the  uniform 
reduction  in  shear  stress  acting  to  produce  seismic  slip  over  a  circular  fault.  For  the  largest  rock  burst  in 
our  data  set,  A a  works  out  to  be  2.7  Mpa. 

2.10  Apparent  Stress. 

The  apparent  stress  aA  is  another  seismic  measure  of  stress  release  (Wyss  and  Brune,  1968).  It  is 
calculated  from  radiated  P  and  S-wave  energy  and  the  seismic  moment: 

aA  =  p  (Ea  +  Ep)  /  M0  (2.10) 

The  apparent  stress  is  proportional  to  the  dynamic  stress  drop,  but  does  not  represent  an  actual  stress 
difference.  For  our  largest  rockburst,  aA  calculates  out  at  317  Pa 

2.11  Dynamic  Stress  drop. 

As  described  above,  the  dynamic  stress  drop  can  be  estimated  from  aRMs,  the  root  mean  square 
acceleration  averaged  over  the  duration  of  the  S  wave  arrival  as  follows: 

°rms  =  3.2  p  R  aRMs  (fs  /  f  max)1/2 

where  fs  =  S  wave  comer  frequency  on  the  displacement  sprectrum,  and 
f»4Ax  =  the  highest  frequency  of  acceleration  reliably  recorded. 

If  the  rupture  process  is  not  complex  and  rupture  velocity  is  constant,  then  Crms  and  Ao  should  be  similar 
in  magnitude.  For  the  largest  rockburst  in  our  study,  Crms  calculates  to  be  l.l  Mpa,  less  than  half  of  the 
static  stress  drop,  Aa,  suggesting  that  asperities  along  the  fault  plane  were  important  in  this  event. 

2.12  Fault  Displacement. 

The  average  slip  across  the  fault  plane  is  useful  for  source  characterization.  For  the  Brune  Model,  the 
average  displacement  u  is  calculated  using  the  definition  of  seismic  moment  rearranged: 


U  =  Mq/  (jlTTo2) 


(2.12) 


For  the  largest  rockburst  in  our  data  set,  the  average  slip  across  the  fault  was  found  to  be  5  mm. 
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'Figure  11.  The  average  day  at  the  Lucky  Friday  Mine.  Although  most  events  occur  at  blasting 
times,  events  can  occur  at  any  time  of  day. 


2.13  Fault  Plane  Solutions  —  Moment  Tensor  Inversion. 

The  largest  rockburst  (Figure  10)  analyzed  in  this  study  yielded  seismograms  with  a  moment  tensor 
solution  involving  a  significant  volumetric  component.  The  moment  tensor  was  decomposed  into 
isotropic  and  deviatoric  components  from  which  the  volumetric  closure  and  total  shear  deformation  were 
estimated  for  comparison  with  in-mine  coseismic  damage.  The  magnitude  3.1  burst  occurred  7.5  hours 
after  routine  daily  blasting  in  the  mine  (Figure  11).  The  event  represented  the  result  of  deep  underground 
mining  combined  with  sliding  on  nearby  bedding  planes  or  faults.  The  implosional  characteristic  may 
prove  to  be  a  useful  parameter  for  fingerprinting  seismic  events  from  deep  hardrock  mines. 

The  moment  tensor  of  the  rock  burst  was  estimated  in  the  far  field  using  a  linear  inversion  method  in  the 
time  domain: 


U  =  G  m 


(2.13) 


where 

U  =  n  sampled  observations  of  ground  displacement  at  the  various  stations, 

G  =  an  n  x  6  matrix  containing  the  Green’s  functions  calculated  by  the  appropriate  algorithm  and 
earth  model,  and 

m  =  (Ml  1,  M12,  M22,  M13,  M23,  M33),  a  vector  containing  the  six  moment  tensor  elements  to  be 
detennined  (Gibowicz  and  Kijko,  1994). 
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Green’s  functions  were  formulated  in  a  simple  manner  using  the  displacement  vectors  for  far  field  P,  SV, 
and  SH  wave  radiation  (Pujol  and  Herrmann ,  1990)  at  each  of  the  five  stations.  A  synchronous  source  in 
the  form  of  a  step  function  was  used  to  estimate  displacements.  Based  on  observations  by  Stiekney  and 
Sprenke  (1993),  we  assumed  a  isotropic  and  homogeneous  velocity  model  with  a  P-wave  velocity  of  5.74 
km/s,  and  S-wave  velocity  of  3.47  km/s.  A  density  of  2.7  was  assumed  for  the  metasedimentary  rocks  in 
the  mine  area. 

The  far-field  contribution  of  a  point  source  to  P  or  S  wave  particle  displacement  at  a  given  location  is 
given  by: 

uc(t)  =  4  ir  pc3  Rc  R'1  MoXt-R/c)  (2.14) 


where  uc  is  displacement  of  the  P,  SV,  or  SH  wave,  Rc  is  the  radiation  pattern,  and  M0’(t)  is  the  time 
history  of  the  seismic  moment  Mo .  Following  McGarr  (1991),  the  time  history  of  each  moment  tensor 
component  is  assumed  to  be: 

Mo*(t)  =  M0  f(l  -cos2  ?rft)  for  0  <t  ^1/f  (2.15) 

0  otherwise 

where  f  is  the  inverse  of  the  displacement  pulse  duration.  Because  we  measure  the  peak  displacement  of 
the  three  waves,  we  use  M0’(t)  =  2  M0  f  at  the  time  of  the  maximum  pulse  displacement  and  Mo'(t)  =  0 
otherwise. 

Figure  12  shows  Greens  functions  for  the  study  area  centered  on  the  Lucky  Friday  mine.  For  the 
simulation,  the  surface  is  considered  flat.  The  nature  of  P,  SV  and  SH  displacement  can  be  readily  seen  in 
comparison  with  the  nine  couples  that  compose  the  seismic  moment  tensor  (Figure  13). 
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Figure  12.  Greens  functions  for  the  study  area.  The  simulated  M3.1  hypocenter  is  at  a  depth  of 

1750  m  in  the  center  of  the  mine  area.  The  accelerometer  stations  are  indicated.  The  red 
areas  are  compressional,  the  blue,  dilatational. 


2.14  Analysis  of  Near-Field  Data. 

For  this  project,  we  investigated  the  seismic  signals  from  a  damaging  rock  burst  in  the  Lucky  Friday 
Mine.  The  local  magnitude  of  the  event  was  3.1  as  measured  at  Butte,  Montana.  It  was  recorded  on  the 
Montana  and  Washington  region  networks.  Predominantly  dilatational  first  motions  for  this  event  from 
our  five  local  stations  and  several  regional  stations  suggested  that  this  event  had  an  implosional 
component.  To  test  this  hypothesis,  we  performed  a  moment  tensor  inversion  using  the  signals  from  five 
strong  motion  accelerometers  that  were  deployed  in  the  immediate  area  of  the  mine  (Figure  9). 

The  seismic  signals  used  in  this  experiment  were  collected  as  part  of  a  local  network  deployment  of  five 
strong  motion  accelerometers  and  one  broadband  seismometer  from  6-1 1-98  to  9-2-98  in  the  immediate 
vicinity  of  the  Lucky  Friday  Mine.  The  damaging  rock  burst  was  the  largest  seismic  event  recorded  in  that 
period.  One  accelerometer  was  located  in  an  abandoned  adit;  the  remaining  stations  were  surface  stations. 
The  data  used  were  the  maximum  displacement  amplitudes  of  the  P,  SV,  and  SH  waves  as  recorded  at  the 
five  three-component  strong  motion  stations  surrounding  the  site  (Figures  14-18).  Linear  trends  in  the 
records  were  identified  and  removed.  Because  of  their  clarity,  the  records  were  not  filtered.  The 
amplitudes  of  the  P,  SH,  and  SV  phases  were  measured  at  each  of  the  five  stations,  and  corrected  for  the 
free  surface  effect.  The  source  location  was  provided  by  the  in-mine  geophone  array  operated  by  the 
mine. 
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Figure  14.  Displacement  amplitudes  at  ATL  for  the  largest  burst  recorded  in  this  study.  The  P  and 
S  picks  are  also  shown.  See  Figure  9  for  locations  of  the  stations. 
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Figure  15.  Displacement  amplitudes  at  DED  for  the  largest  burst  recorded  in  this  study.  The  P 
and  S  picks  are  also  shown.  See  Figure  9  for  locations  of  the  stations. 
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Figure  16.  Displacement  amplitudes  at  MOR  for  the  largest  burst  recorded  in  this  study.  The  P 
and  S  picks  are  also  shown.  See  Figure  9  for  locations  of  the  stations. 
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Figure  17.  Displacement  amplitudes  at  GOL  for  the  largest  burst  recorded  in  this  study.  The  P  and 
S  picks  are  also  shown.  See  Figure  9  for  locations  of  the  stations. 
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Displacement  meters 


To  estimate  the  moment  tensor  elements,  we  followed  a  least  squares  procedure  outlined  by  MeGarr 
(1992),  We  found  the  following  moment  tensor  solution: 


M 


-0.0743  -0.1041 

-0.1373 


-0.3102 

-0.0738 

-1.0000 


x  5.03  x  10  13  N-m 


(2.16) 


The  plots  of  surface  displacements  for  this  solution  in  the  mine  area  show  that  the  mechanism  is  an 
excellent  example  of  a  shear-implosional  failure.  The  P  displacements  (Figure  19)  are  dominated  by  the 
large  negative  M33  tensor  component.  This  results  in  dilatational  first-motions  on  all  seismic  stations. 
The  implosional  nature  of  the  rock  burst  is  clearly  evident  on  this  figure.  The  stereographic  projection  of 
the  P  first  motions  (Figure  20)  shows  all  dilations  except  for  a  small  region  of  compressive  first  motions, 
The  SV  and  SH  displacements  (Figures  21  and  22)  appear  to  be  dominated  by  a  negative  M13 
component,  suggesting  slip  on  faults  that  cut  through  the  mine  area.  The  shear  nature  of  the  event  is  very 
evident  on  the  S  displacements. 


Figure  19.  The  P  displacements  for  the  largest  burst  observed  the  shown  on  a  plan  view  of  the  mine 
area.  The  implosional  (negative  M33)  nature  of  the  event  is  clearly  evident. 
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Figure  20.  The  P  displacement  polarity  plotted  on  the  lower  hemisphere  of  an  equal-area 
stereographic  projection. 
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Figure  21.  The  SV  displacements  for  the  largest  burst  observed  shown  on  a  plan  view  of  the  mine 
area.  The  shear  (negative  M13)  nature  of  the  event  is  evident. 
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Figure  22.  The  SH  displacements  for  the  largest  burst  observed  shown  on  a  plan  view  of  the  mine 
area.  The  shear  (negative  M13)  nature  of  the  event  is  apparent. 
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To  further  interpret  the  moment  tensor  elements,  we  first  followed  the  procedure  outlined  by  McGarr 
(1992)  in  his  study  of  South  African  rock  bursts.  The  amount  of  coseismic  closure  (AV)  due  to  the 
implosion  can  be  computed  using  the  trace  of  the  moment  tensor  along  with  the  estimated  elastic 
properties  of  the  metasedimentary  quartzites  in  the  mine  area: 

AV  =  tr(Mij)/(3X+2p)  =  (-6.09  x  10  13  N-m)/(1.63  x  10  11  Pa)  (2.17) 

r 

or  =  -379  m3  of  coseismic  volume  change. 

To  compute  the  shear  deformation,  the  moment  tensor  is  first  diagonalized:  - 


M(diag)  = 


0.0536  0  0 

0  -0.1603  0 

0  0  -1.1049 


x  5.03  x  10  13  N-m 


(2.18) 


The  Eigenvectors,  representing  directions  of  the  seismic  P,  I,  and  T  axes,  respectively  are: 


(  0.8887  0.3499  -0.2963) 
(-0.3942  0.9131  -0.1043) 
(-0.2340  -0.2095  -0.9494) 


(2.19) 


The  deviatoric  component  is  estimated  by  subtracting  the  volumetric  component: 


M(dev)  =  M(diag)  - 1/3  tr(Mij) 


(2.20) 


M(dev)  = 


0.4566  0  0 

0  0.2427  0 

0  0  -0.7019 


x  5.03  x  10  13  N-m 


Decomposing  into  a  major  and  minor  double-couple: 


0.70 

-0.24 

M(dev)  = 

0 

-0.70 

+ 

0.24 

0 

(2.21) 


For  the  major  double-couple,  the  P  axis  is  vertical  and  the  T  axes  is  east .  For  the  minor  double  couple, 
the  P  axis  is  east  and  the  T  axis  is  north..  These  agree  well  with  the  dominant  east-northeasterly  direction 
of  tectonic  extension  in  the  area. 


* 
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The  total  shear  deformation: 


2  AD  =  {Ml  (dev)  +  M2(dev)}/p 

or  2  AD  =  [(0.46+0.24)  x  5.03  x  10  13  N-m]  /  (3.76  x  1010  Pa) 

2  AD  =  909  m3 , 2.4  times  larger  than  the  volumetric  reduction. 


(2.22) 
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An  alternative  approach  to  interpret  the  moment  tensor  solution  is  to  apply  a  crack  closure  (CC) 
mechanism.  This  is  the  negative  of  the  tension  crack  (TC)  suggested  by  Aki  and  Richards  (1980): 


M(crack)  = 


vx 

0 

0 


0  0 

vx  0 

0  x 


(2.23) 


where  v-A./(X  +  2p)  ~  1/3,  x  -  (X  +  2  p)  AV,  and  M(crack)  is  the  moment  tensor  of  a  crack  with 
change  in  volume  AV  =  u  A,  where  u  is  the  crack  closure  displacement  and  A  is  the  area  of  the  crack. 

To  implement  this  interpretation,  we  begin  with  the  diagonalized  moment  tensor: 


M(diag)  = 


Ml  0  0 

0  M2  0 

0  0  M3 


(2.24) 


We  attempt  a  decomposition  as  follows  to  break  the  diagonalized  moment  tensor  into  a  crack  closure,  a 
double  dipole  M(dd),  and  an  error  term  M(e). 


4 


V 


M(diag)=  M(crack)  +  M(dd)  +  M(e) 

M(dd)  can  take  three  forms,  depending  on  the  orientation  of  the  double  dipole: 


M’ 

0 

0 

0 

0 

0 

-M’ 

0 

0 

0 

-M’ 

0 

or 

0 

M’ 

0 

or 

0 

0 

0 

0 

0 

0 

0 

0 

-M’ 

0 

0 

M’ 

(2.25) 


(2.26) 


For  our  data,  the  third  form  above  gave  the  best-fit  result.  We  obtain  three  simultaneous  equations. 

Ml  =  x/3  -  M’ 

M2  =  x/3  (2.27) 

M3  =  x  +  M’ 
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The  least  square  solution  to  these  equations  yields  an  estimate  of  x  and  M’  and  the  error  M(e). 
For  our  data,  we  obtain  x  =  -0.7546  x  5.03  x  1013  N-m  and  M’  =  0.3276  x  5.03  x  1013  N-m. 


M(crack)  = 


-0.2513  0  0 

0  -0.2515  0 


x  5.03  xlO'3  N-m 


0  0  -0.7546 


(2.28) 


M(dc)  = 


0.3276  0 

0  0 

0  0 


0 

0 

-0.3276 


x  5,03  xlO13  N-m 


(2.29) 


M(«) 


0.0227  0  0 

0  -0.0912  0 

0  0  0.0227 


x  5.03  x  1013  N-m 


(2.30) 
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This  result,  with  an  overall  RMS  error  of  0.0557  x  5.03  x  1013  N-m,  is  the  best  fit  of  the  three  cases 
considered  for  the  double  couple.  The  other  two  cases  produced  errors  4  to  5  times  larger.  Comparing 
with  the  Eigenvectors  above,  this  double  couple  has  an  east-northeast  T-axis  and  a  near  vertical  P-axis. 
This  solution  is  consistent  with  the  regional  stress  direction,  and  it,  like  the  McGarr  model,  produces  the 
anticipated  normal  slip  predicted  by  the  negative  M13  component  of  the  stress  tensor  solution. 

The  volume  reduction  in  the  mine  AV  —  x/(  X  +  2  \i)  works  out  to  be  -420  m3 .  The  shear  deformation  due 
to  the  double  couple  E  AD  =  M7/r  is  as: 

E  AD  =  [(0.3276)  x  5.03  x  10  13  N-m]  I  (3.76  x  10'°  Pa)  (2.31) 

£  AD  =  350  m3 ,  or  about  83%  of  the  volumetric  reduction. 

It  is  interesting  to  note  that  the  M(e)  error  term  from  the  least  squares  fit  is  similar  in  form  to  a 
compensated  linear  vector  dipole  (CVLD): 


M(CVLD)  = 


-x72  0  0 

0  -x’  0 

0  0  -x72 


(2.32) 


except  that  our  result  appears  to  have  one  dipole  of  strength  4  rather  than  strength  2  relative  to  the  others. 
At  any  rate  the  M(e)  tensor  is  negligible  in  our  result. 


* 
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2.15  A  Template  for  Shear-Implosional  Rockbursts. 

The  moment  tensor  estimates  for  the  rock  burst  analyzed  here  were  found  to  be  dominated  by  negative 
diagonal  components,  indicating  a  dominantly  implosional  source.  However,  the  source  is  not  completely 
isotropic.  Our  results  are  similar  to  those  obtained  for  a  group  of  events  studied  by  McGarr  (1992)  from 
deep  hard  rock  mining  districts  in  South  Africa,  suggesting  that  implosional  events  may  be  characteristic 
of  a  type  of  event  that  may  fingerprint  such  mines  for  CTBT  purposes.  Further  work  needs  to  be  done  to 
'  evaluate  the  effect  of  such  implosional  events  at  distances  useful  for  treaty  verification  monitoring. 

Regional  waveforms  were  obtained  for  nearby  earthquakes  and  explosions  that  have  occurred  in  the 
<  region  (Figure  23)  and  for  the  larger  events  analyzed  in  this  study  (Figure  24).  These  data  were  used  in 

the  second  phase  of  this  contract  to  evaluate  the  ability  to  discriminate  between  such  events  and  shear-  or 
implosion-dominated  mine-induced  events  using  regional  waveforms. 


Northern  Idaho  Mine  Tremor  and  Earthquake  at  WALA 


Figure  23.Waveforms  at  CSN.  seismic  station  VVALA  from  a  3.6  mine  tremor  on  Feb.  9, 1996 
(upper)  and  from  a  recent  northern  Idaho  earthquake  (lower).  Approximate  group 
velocities  of  7.8, 6.1,  and  3.7  km/s  are  shown. 
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The  key  to  fingerprinting  many  of  the  rock  bursts  from  deep  underground  mines  such  as  those  in  the 
Silver  Valley  is  to  make  use  of  the  implosional  component  of  the  bursts  that  seems  to  be  associated  with 
the  majority  of  the  bursts  and  the  consequent  lack  of  Love  waves  associated  with  such  sources.  For 
example,  Pechman  et  al.  (1995)  have  shown  that  five  broadband  stations  at  distances  from  300  and  1000 
km  were  sufficient  to  resolve  the  focal  mechanism  and  seismic  moment  of  the  implosional  ML5.1  Trona 
Mine  event.  In  the  present  study,  we  have  shown  that  the  source  characteristics  of  the  rockbursts  in  the 
Coeur  d'Alene  Mining  District  are  similar  to  those  of  the  Trona  Mine  event.  Thus  the  methods  used  to 
characterize  the  Trona  Mine  event  are  also  applicable  to  deeper  underground  events  as  well. 

The  Trona  Mine  event  was  characterized  by  modeling  its  long  period  waveforms  using  a  regional 
velocity  model  in  conjunction  with  methodology  as  described  by  Walter  (1993).  The  seismograms 
evaluated  were  instrument-corrected  displacement  records  band  passed  filtered  from  20  to  40  seconds. 
The  seismograms  showed  little  long  period  Pnl  or  Love  wave  amplitude.  Pnl  and  Love  wave  energy  is 
characteristic  of  earthquakes,  but  not  explosions  or  mine  collapses.  Sufficient  broadband  stations  are 
necessary,  however,  to  ensure  that  azimuthal  coverage  is  sufficient  to  avoid  Love  Wave  nodes  at 
particular  azimuths. 

The  horizontal  tension  crack  model  appears  to  at  least  partially  characterize  both  shallow  and  deep  mine 
collapses.  The  collapsing  crack  model  produces  all  dilatational  first  arrivals,  and  for  the  case  of  the  Trona 
Mine  events,  fits  the  waveforms  at  distant  stations  better  than  a  pure  implosional  model  or  an  earthquake 
model.  For  the  shear  implosional  mechanism  apparent  in  the  Coeur  d'Alene  District  events,  a  combination 
of  the  collapsing  crack  model  and  a  double  couple  needs  to  be  employed  to  better  fit  the  waveforms  at 
regional  distances.  Figure  24  compares  a  rock  burst  and  earthquake  record  from  north  Idaho.  Both 
events  clearly  have  well  developed  Love  waves,  so  a  simple  implosional  or  collapsing  crack  model  will 
fail  to  adequately  model  the  waveforms.  A  shear  implosional  source  would  have  to  be  employed  instead. 

Seismic  events  related  to  mining  are  commonly  thought  to  be  of  two  distinct  types.  Type  1  events  are 
caused  by  implosion;  type  2  are  tectonic  events  involving  fault  slip  triggered  somehow  by  mining.  The 
deep  underground  rock  bursts  in  the  Silver  valley  appear  to  be  a  hydrid  event  involving  both  implosion 
and  shear  along  failure  planes. 
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